function DataPlot(fileNE)%filenamec0,fresult1,filenamec1,fresult2,filenamec2
%输入文件名，相应的设置再下方的备注中有

data_matrix = readmatrix(fileNE);
N=data_matrix(:,1);
E=data_matrix(:,2);
plot(N,E);
title('有零速修正小推车平面轨迹');
xlabel('北方向(m)');
ylabel('东方向(m)');



%{
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%示例数据误差
data_matrix = readmatrix(filenamec0);
delta_phi=data_matrix(:,2);
delta_lam=data_matrix(:,3);
delta_h=data_matrix(:,4);
delta_vn=data_matrix(:,5);
delta_ve=data_matrix(:,6);
delta_vd=data_matrix(:,7);
delta_Phi=data_matrix(:,8);
delta_theta=data_matrix(:,9);
delta_psi=data_matrix(:,10);
for i=1:length(delta_psi)
    if(delta_psi(i)<-180)
        delta_psi(i)=delta_psi(i)+360;
    end
end
x=(1:length(delta_phi))';

figure;
subplot(3,1,1); %位置
plot(x,delta_phi);
title('示例数据纬度误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2); 
plot(x,delta_lam);
title('示例数据经度误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,delta_h);
title('示例数据高程误差');
xlabel('历元数');
ylabel('m');
figure;
subplot(3,1,1); %速度
plot(x,delta_vn);
title('示例数据北向速度误差');
xlabel('历元数');
ylabel('m/s');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,delta_ve);
title('示例数据东向速度误差');
xlabel('历元数');
ylabel('m/s');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,delta_vd);
title('示例数据地向速度误差');
xlabel('历元数');
ylabel('m/s');
figure;
subplot(3,1,1); %姿态
plot(x,delta_Phi);
title('示例数据横滚角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,delta_theta);
title('示例数据俯仰角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,delta_psi);
title('示例数据航向角误差');
xlabel('历元数');
ylabel('deg');

%}
%%
%{
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%无零速修正解算结果
data_matrix = readmatrix(fresult1);
phi=data_matrix(:,2);
lam=data_matrix(:,3);
h=data_matrix(:,4);
vn=data_matrix(:,5);
ve=data_matrix(:,6);
vd=data_matrix(:,7);
Phi=data_matrix(:,8);
theta=data_matrix(:,9);
psi=data_matrix(:,10);
x=(1:length(vn))';

figure;
subplot(3,1,1); %位置
plot(x,phi);
title('无零速修正纬度变化图');
xlabel('历元数');
ylabel('经纬度(deg)');
subplot(3,1,2); 
plot(x,lam);
title('无零速修正经度变化图');
xlabel('历元数');
ylabel('经纬度(deg)');
subplot(3,1,3); 
plot(x,h);
title('无零速修正高程变化图');
xlabel('历元数');
ylabel('高程(m)');
figure;
subplot(3,1,1); %速度
plot(x,vn);
title('无零速修正北向速度变化图');
xlabel('历元数');
ylabel('北向速度(m/s)');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,ve);
title('无零速修正东向速度变化图');
xlabel('历元数');
ylabel('东向速度(m/s)');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,vd);
title('无零速修正地向速度变化图');
xlabel('历元数');
ylabel('地向速度(m/s)');
figure;
subplot(3,1,1); %姿态
plot(x,Phi);
title('无零速修正横滚角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,theta);
title('无零速修正俯仰角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,psi);
title('无零速修正航向角变化图');
xlabel('历元数');
ylabel('deg');
%}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%不零速修正解算结果误差
%{
data_matrix = readmatrix(filenamec1);
delta_phi=data_matrix(:,2);
delta_lam=data_matrix(:,3);
delta_h=data_matrix(:,4);
delta_vn=data_matrix(:,5);
delta_ve=data_matrix(:,6);
delta_vd=data_matrix(:,7);
delta_Phi=data_matrix(:,8);
delta_theta=data_matrix(:,9);
delta_psi=data_matrix(:,10);
for i=1:length(delta_psi)
    if(delta_psi(i)<-180)
        delta_psi(i)=delta_psi(i)+360;
    end
end
x=(1:length(delta_phi))';

figure;
subplot(3,1,1); %位置
plot(x,delta_phi);
title('无零速修正纬度误差');
xlabel('历元数');
ylabel('经纬度(deg)');
subplot(3,1,2); 
plot(x,delta_lam);
title('无零速修正经度误差');
xlabel('历元数');
ylabel('经纬度(deg)');
subplot(3,1,3); 
plot(x,delta_h);
title('无零速修正高程误差');
xlabel('历元数');
ylabel('高程(m)');
figure;
subplot(3,1,1); %速度
plot(x,delta_vn);
title('无零速修正北向速度误差');
xlabel('历元数');
ylabel('北向速度(m/s)');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,delta_ve);
title('无零速修正东向速度误差');
xlabel('历元数');
ylabel('东向速度(m/s)');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,delta_vd);
title('无零速修正地向速度误差');
xlabel('历元数');
ylabel('地向速度(m/s)');
figure;
subplot(3,1,1); %姿态
plot(x,delta_Phi);
title('无零速修正横滚角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,delta_theta);
title('无零速修正俯仰角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,delta_psi);
title('无零速修正航向角误差');
xlabel('历元数');
ylabel('deg');

%}
%%
%{
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%零速修正解算结果
data_matrix = readmatrix(fresult2);
phi=data_matrix(:,2);
lam=data_matrix(:,3);
h=data_matrix(:,4);
vn=data_matrix(:,5);
ve=data_matrix(:,6);
vd=data_matrix(:,7);
Phi=data_matrix(:,8);
theta=data_matrix(:,9);
psi=data_matrix(:,10);
x=(1:length(vn))';

figure;
subplot(3,1,1); %位置
plot(x,phi);
title('零速修正纬度变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2); 
plot(x,lam);
title('零速修正经度变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,h);
title('零速修正高程变化图');
xlabel('历元数');
ylabel('高程(m)');
figure;
subplot(3,1,1); %速度
plot(x,vn);
title('零速修正北向速度变化图');
xlabel('历元数');
ylabel('北向速度(m/s)');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,ve);
title('零速修正东向速度变化图');
xlabel('历元数');
ylabel('东向速度(m/s)');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,vd);
title('零速修正地向速度变化图');
xlabel('历元数');
ylabel('地向速度(m/s)');
figure;
subplot(3,1,1); %姿态
plot(x,Phi);
title('零速修正横滚角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,theta);
title('零速修正俯仰角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,psi);
title('零速修正航向角变化图');
xlabel('历元数');
ylabel('deg');
%}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%零速修正解算结果误差
%{
data_matrix = readmatrix(filenamec2);
delta_phi=data_matrix(:,2);
delta_lam=data_matrix(:,3);
delta_h=data_matrix(:,4);
delta_vn=data_matrix(:,5);
delta_ve=data_matrix(:,6);
delta_vd=data_matrix(:,7);
delta_Phi=data_matrix(:,8);
delta_theta=data_matrix(:,9);
delta_psi=data_matrix(:,10);
for i=1:length(delta_psi)
    if(delta_psi(i)<-180)
        delta_psi(i)=delta_psi(i)+360;
    end
end
x=(1:length(delta_phi))';

figure;
subplot(3,1,1); %位置
plot(x,delta_phi);
title('零速修正纬度误差');
xlabel('历元数');
ylabel('纬度(deg)');
subplot(3,1,2); 
plot(x,delta_lam);
title('零速修正经度误差');
xlabel('历元数');
ylabel('经度(deg)');
subplot(3,1,3); 
plot(x,delta_h);
title('零速修正高程误差');
xlabel('历元数');
ylabel('高程(m)');
figure;
subplot(3,1,1); %速度
plot(x,delta_vn);
title('零速修正北向速度误差');
xlabel('历元数');
ylabel('北向速度(m/s)');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,delta_ve);
title('零速修正东向速度误差');
xlabel('历元数');
ylabel('东向速度(m/s)');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,delta_vd);
title('零速修正地向速度误差');
xlabel('历元数');
ylabel('地向速度(m/s)');
figure;
subplot(3,1,1); %姿态
plot(x,delta_Phi);
title('零速修正横滚角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,delta_theta);
title('零速修正俯仰角误差');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,delta_psi);
title('零速修正航向角误差');
xlabel('历元数');
ylabel('deg');
%}
%%
%{
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%参考值画图
data_matrix = readmatrix(fTrue);
phi=data_matrix(:,3);
lam=data_matrix(:,4);
h=data_matrix(:,5);
vn=data_matrix(:,6);
ve=data_matrix(:,7);
vd=data_matrix(:,8);
Phi=data_matrix(:,9);
theta=data_matrix(:,10);
psi=data_matrix(:,11);
x=(1:length(vn))';

figure;
subplot(3,1,1); %位置
plot(x,phi);
title('参考纬度变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2); 
plot(x,lam);
title('参考经度变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,h);
title('参考高程变化图');
xlabel('历元数');
ylabel('高程(m)');
figure;
subplot(3,1,1); %速度
plot(x,vn);
title('参考北向速度变化图');
xlabel('历元数');
ylabel('m/s');
subplot(3,1,2); % 2行1列的第1个位置
plot(x,ve);
title('参考东向速度变化图');
xlabel('历元数');
ylabel('m/s');
subplot(3,1,3); % 2行1列的第1个位置
plot(x,vd);
title('参考地向速度变化图');
xlabel('历元数');
ylabel('m/s');
figure;
subplot(3,1,1); %姿态
plot(x,Phi);
title('参考横滚角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,2);
plot(x,theta);
title('参考俯仰角变化图');
xlabel('历元数');
ylabel('deg');
subplot(3,1,3); 
plot(x,psi);
title('参考航向角变化图');
xlabel('历元数');
ylabel('deg');
%}

end